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ABSTRACT 

In the theory of star formation, the first hydrostatic core (PHSC) phase is a critical step in which a condensed object emerges from a 
prestellar core. This step lasts about one thousand years, a very short time compared with the lifetime of prestellar cores, and therefore 
is hard to detect unambiguously. We present IRAM Plateau de Bure observations of the Barnard lb dense molecular core, combining 
detections of H 2 CO and CH 3 OH spectral lines and dust continuum at 2.3” resolution (~ 500 AU). The two compact cores Blb-N and 
B Ib-S are detected in the dust continuum at 2mm, with fluxes that agree with their spectral energy distribution. Molecular outflows 
associated with both cores are detected. They are inclined relative to the direction of the magnetic field, in agreement with predictions 
of collapse in turbulent and magnetized gas with a ratio of mass to magnetic flux somewhat higher than the critical value, p ~ 2-1. 
The outflow associated with Blb-S presents sharp spatial structures, with ejection velocities of up to 7 kms“* from the mean 
velocity. Its dynamical age is estimated to be ~ 2000 yr. The Blb-N outflow is smaller and slower, with a short dynamical age of 
~ 1000 yr. The Blb-N outflow mass, mass-loss rate, and mechanical luminosity agree well with theoretical predictions of PHSC. 
These observations confirm the early evolutionary stage of Blb-N and the slightly more evolved stage of Blb-S. 
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00 1- Introduction 

In the current theory of low-mass star formation, the first con- 
'sD densed object to be formed in the collapse of the prestellar 
cores is the first hydrostatic core (FHSC). During this rather 
short-lived phase of less than about a thousand years, the col- 
lapse slows down until the core temperature rises above about 
2000 K, when molecular hydrogen starts to dissociate dLarsonl 
T—I Il969ll . Predictions of the spectral energy distribution of FH- 
rG SCs as well as their density structure and kinematics are avail- 
. able thanks to magnetohyd rodynamic (MHD) simulations (e.g. 
ICommercon et al.|l2012allbh . Their observational characteristics 
$_i are intermediate between those of a prestellar core and a class 
0 protostar: i) a spectral energy distribution (SED) peaking be¬ 
yond ~ 100 pm similar to that of a prestellar core, and ii) 
a compact structure of ~ 500 AU. These features are found 
in very low luminosity objects (VeLLOs), which are compact 
sour ces characterized by an internal luminosity lower than 0.1 
Lq ddi Francesco et all l2007l) . They can be class 0 or class I 
young stellar objects (YSOs), proto-brown dwarfs, or FHSC can¬ 
didates. VeLLOs exhibit contrasting outflow properties, from ex¬ 
tended and collimated lobes dow n to faint, slow, and compact 
flows {e.g. iDunham et al.l d201 ih l. The characterization of an 


* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (Prance), MPG 
(Germany) and IGN (Spain). 


FHSC candidate requires studying its SED and the gas kinemat¬ 
ics of its envelope and molecular outflow. Another important the¬ 
oretical prediction of the FHSC stage is the presence of a com¬ 
pact, slow, and poorly collimated outflow, whose properties de¬ 
pend on the strength and orientation of t he magnetic field relative 


to the rotation axis of the prestellar core (iHennebelle & Fromang 

2008 

: iTomida et al.l 120101; ICiardi & Hennebellel 20101; Machida 

2014 

). 


As of today, several FHSC candidates have been found 
in star-forming regions in the Perseus molecular cloud. Of 
these, the Barnard 1 region has attracted attention because of 
it has dense cores at different evolutionary stages: Bla and 
Blc a re each hosting a class 0 YSO with developed out¬ 
flows (iHatchell et aLlUOPTl) . while Bib has recently been pro¬ 
posed to host two FHS Cs based on the SED analysis of Her- 
schel and Spitzer data (iPezzuto et al.l l2012ll . In detail, three 
remarkable sources are detected in Bib: an infrared source 
det ected by Spitzer, with strong absorptio n from ices Blb- 
W (l.lOrgensen et al.ll2006l : lEvans et '^120091) . and two compact 
(sub)millimeter continuum sources Blb-N and Blb-S that are 
the FHSC candidates. The three sources are deeply embedded in 
the surrounding protostellar core, which seems essentially unaf¬ 
fected by them, as shown by its high column density, AfH?) ~ 
10^^ cm~^ and low kinetic temperature, = 12 K (iLis et al.l 

1201 nil . 
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Because of the nearby Bla and Blc YSOs, it has been diffi- 
cu lt to associate molecu lar outflows with either Blb-N or Blb- 
S. iHatchell et al.l (l2007h found high-velocity emission close to 
Blb-S in their CO(3-2) survey of the region, but without a 
clear association with either source. Broad line profiles and line 
wings have also be en detected in some species such a s CH^OH 
( Oberg et al.l2010l) and the excited transition NH 3 ( 3 , 3 ) dLis et al.l 
l2010l) ywithout assignment to any of the Bib sources. New SMA 
observation of CO(2-l) at a resolution of 6.2" x 3.9" unambigu¬ 
ously detected mole cular outflows with the Blb-S and Blb-W 
and possibly Blb-N (iHirano & Liull20f^ . 

To better understand the nature of the three objects em¬ 
bedded in Bib, we have obtained 2 " resolution data with the 
IRAM Plateau de Bure interferometer, targeting methanol and 
formaldehyde lines in the 2 mm spectral window. Following 
iHirano & Liul (|2014|) . we adopt a distance of 230 pc for Barnard 
1, implying that 2" corresponds to ~ 500 AU. The observations 
are presented in Sect.|2]and are discussed in Sect.[3]with partic¬ 
ular emphasis on the molecular outflows. 


2. Observations 

Table 1 summarizes the interferometric and single-dish observa¬ 
tions. The calibration and the joint imaging and deconvolution 
processes are described in Appendix lAl 

2.1. Interferometric observations 

For the interferometric observations, the lower sideband of the 
2 mm receivers was tuned at 145.3 GHz. The WIDEX backend 
yielded a total bandwith of 4 GHz per polarization at a spec¬ 
tral channel spacing of 2 MHz. The intermediate frequency was 
further split into two 1 GHz chunks centered around 144.5 and 
145.3 GHz. The eight windows (40MHz) of the high spectral 
resolution correlator were then centered around potential lines 
inside the previously defined WIDEX chunks. This yielded spec¬ 
tra with a 78 kHz channel spacing that we further smoothed to 
reach a spectral resolution of 0.2 km s ' (or typically 97 kHz). 
We targeted H 2 CO( 2 o ,2 - lo.i) at 145.60295 GHz and the 
band of CH 30 H (3 - 2 ) lines with frequencies ranging from 
145.0937 GHz to 145.1318 GHz. We observed a mosaic of seven 
pointings that followed a hexagonal compact pattern with near¬ 
est neighbors separated by half the primary beam, providing a 
roughly circular held of view of 70" diameter at 2 mm. 

2.2. Single-dish observations 

The IRAM-30m data were simultaneously observed at 3 and 
2 mm with a combination of the EMIR receivers and the Fourier 
transform spectrometers, which yields a bandwidth of 1.8 GHz 
per sideband per polarization at a channel spacing of 49 kHz. The 
mixers were tuned at 85.55 and 144.90 GHz. We used the on-the- 
fly scanning strategy with a dump time of 0.5 s and a scanning 
speed of 10.9"/s to ensure a sampling of at least three dumps per 
beam at 17.8", the resolution at 145 GHz. The 175" x 175" map 
was covered using successive orthogonal scans along the RA 
and DEC axes. The separation between two successive rasters 
was 6.5" (~ AjlD) to ensure Nyquist sampling at the highest 
observed frequency. A common reference position located at off¬ 
sets (200", 150") was observed for 10 s every 77.5 s. The typical 
IRAM-30m position accuracy is ~ 3". 


Table 2. Physical parameters derived from the CH 3 OH emission 


Position" 

J Tdv'’ 

Kkms^' 

Ne^ 

10 '» cm -2 

KH 2 ) 

10 ^ cm“^ 

Tk 

K 

[CH 3 OH] 

lo-x 

Blb-S 

( 6 ,- 2 ) 

17.9 

48 

6 

30 

20 

(7,-2) 

12.4 

50 

9 

70 

13 

(-9,6) 

1.9 

6 

6 

30 

2.5 

(-8,5) 

5.9 

13 

4 

50 

8 

(- 8 , 10 ) 

2.9 

7.2 

3 

30 

6 

(-7,1) 

4.9 

14 

6 

30 

5 

Blb-N 

(-5,1) 

0.9 

2 

6 

20 

0.8 


Notes. Offsets are given in arcsec relative to each young stellar ob¬ 
ject. Integrated intensity of the E-CHsOHfS-i - 2_i) line. Ne = 
A(E-CH30H); AfCHjOH) = A(E-CH30H) -r A(A-CH30H) ~2xNe. 


Table 3. Physical properties of the molecular outflows 


Parameter 

Unit 

BIb-N 

Blue" Red* 

Blb-S 

Blue" Red* 

Max. velocity 

km s~‘- 

4.5 

4 

7.6 

6.6 

Size 

10 » AU 

1.3 

0.6 

3 

3.3 

Dyn. time 

10» yr 

2.0 

1.0 

1.9 

2.5 

Mass 

10-3 Mo 

1.2 

0.4 

3.0 

1.9 

Momemtum 

10“^ Mq km s“> 

3.6 

1.3 

13 

6.4 

Mass-loss rate 

IQ-^Moyr-* 

5.3 

3.8 

16 

7.5 

Mom. flux 

10“^ Mq km s“> yr“' 

1.5 

1.0 

6.7 

2.5 

Mech. luminosity 

10-3 Lo 

0.8 

0.5 

8.1 

3.1 

Notes. <"> 0.5 - 5.5 

kms ', ® 7.5 - 12.5 kms ' 

(c) .4 

- 5.5 kms 

-1 


<''>7.5-15 kms-'. 



Fig. 2. Examples of CH3OH and H2CO spectra towards the Blb-S and 
Blb-N outflows. The H 2 CO( 2 o 2 - loi) line is shown in black. The four 
detected CH3OH lines are shown in color: blue 3o - 2q A, cyan 3_i - 2_i 
E, green 3o-2o E, and purple 32-22 E. They have been vertically shifted 
by 5 K for clarity. 

3. Results 

The median noise level achieved over the mosaic is 0.24 K (Tmb) 
in 0.2 km s ' channels and for a typical resolution of 2.3". 

3.1. Geometrical properties of the outflows 

Eigure [1] shows the H 2 CO( 2 o ,2 - lo.i) and CH 3 OH integrated 
intensity maps and the mean velocity held around Blb-N and 
Blb-S. The line profiles at selected positions are shown in 
Eig. 13 Blb-W is associated with a secondary peak of the 
H 2 CO and CH 3 OH emission, with narrow line profiles. We did 
not detect any high-velocity emission that might be associated 
with Blb-W or with the high-velocity CO knot detected by 
IHirano & Liul (l2014h . The region between Blb-N and Blb-S 
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Fig. 1. Left: Integrated intensity (from 0 to 14 kms^') and velocity field of the H 2 CO( 2 o ,2 - lo,i) line. Purple contours show the continuum 
emission at 145 GHz. The plus signs indicate the positions of B Ib-N and B Ib-S. The star shows the position of BIb-W. The black ellipse shows 
the beam size. The small panels show the centroid velocity towards Blb-N (top) and Blb-S (bottom) in kms^'. Right: Integrated intensity images 
(from 1.5 to 12.5 kms^*) of the four CH 3 OH lines detected towards Bib: from left to right and top to bottom: 3o - 2o A, 3_i - 2_i E, 3o - 2o E 
and 32 - 1i E. The intensity scale in Kkms^' is given by the color bar on the right. Contours of the 3o - 2o A line are superimposed in all plots, 
ranging from 1 to 20 K km s“ * with a step of 3 K km s“ *. Offsets are given relative to the pointing center given in Table 1. 

Table 1. Observation parameters. The projection center for all observations is 0-2000 = 03''33"'21.300% ^2000 = 31°07'34.00”, a position in between 
Blb-N and Blb-S. 


Molecule 

Frequency 

GHz 

Instr. 

Config. 

Beam 

// 

PA 

Vel. res. 
km s“' 

Int. Time" 
hr 

Tsys 

K 

Noise" 

mK 

Obs. date 

CH 30 H‘' 

145.1 

PdBI 

C&D 

2.30 X 2.24 

128 

0.16 

8.8/29 

160 

313 

Aug. & Oct. 2013 

H 2 C 0 ‘' 

145.6 

PdBI 

C&D 

2.34 X 2.24 

108 

0.20 

8.8/29 

160 

288 

Aug. & Oct. 2013 


Molecule 

Frequency 

GHz 

Instr. 

Feff 

Beff 

Res. 

Res. 

kms“* 

Int. Time" 
hr 

Tsys 

K 

Noise 

mK 

Obs. date 

CH 3 OH 

145.1 

E150 

0.93 

0.74 

17.9 

0.16 

14 

99 

44 

Dec. 3rd & 4th, 2013 

H 2 CO 

145.6 

E150 

0.93 

0.74 

17.8 

0.20 

14 

99 

40 

Dec. 3rd & 4th, 2013 


Notes. Listed as on-source time/telescope time. Evaluated at the mosaic phase center (the noise steeply increases at the mosaic edges after 
correction for primary beam attenuation). Telescope time. FITS files of the H 2 CO and CH 3 OH mosaics are available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 


presents lower CH 3 OH em ission despite the high dust column 
density (iDaniel et al.ll2013l) . a behavior characteristic of deple¬ 
tions due to freezing of molecules onto grains, and supported by 
the detection of high column densitie s of ices including CH 3 OH 
towards Blb-W dBoogert et al.ll200^ . 

The main result is the clear detection of molecular outflows 
in both H2CO and CH3OH around Blb-N and Blb-S. The red 
and blue lobes of the Blb-N outflow are marginally resolved, 
but the excellent correspondence of the H2CO and CH3OH line 
profiles (Fig. EJ, and the alignment with the Blb-N position sup¬ 
port the assignment. The outflow of Blb-S is more extended, 
with its red and blue lobes detected over ~ 10". They are well 
separated, indicating that the inclination angle along the line of 
sight, i, is located between a and 90° - a, where a is the open¬ 
ing angle of either lobe. Using the more extended red lobe, we 
determine a = 23° + 8°, leading to 17° < i < 73°. The red and 
blue lobes of the molecular outflows associated with Blb-N and 
Blb-S appear on different sides of these sources, indicating that 
the outflows have different spatial orientations. 

The line-of-sight component of the magnetic field towards 
Barnard 1 was determined to be B// = -27 + 4 /uG from OH 
Zeeman observations in a 2.9’ beam, toward a position close 
to IRAS 03301-H3057 and ~ 43" from Bib dCioodman et all 


Il989h . The projected direction of the magnetic field is known 
from polarimetric observations of the 850 p. m dust thermal emis¬ 
sion \vith the JCMT at 1 5" resolution dMatthews & Wils^ 
I 2 OO 2 I : [Matthews et alJl200^ . It presents a regular pattern with 
a mean N-S orientation at PA = 1° + 19°. Hence we conclude 
that the magnetic field lies in the vertical plane containing the 
line of sight and the declination axis, its exact orientation de¬ 
pending on the relative strengths of its parallel and perpendicu¬ 
lar components. Therefore, the outflows associated with Blb-N 
and Blb-S are not aligned with the magnetic field, but present 
a significant angle with its mean direction. Using the available 
constraints on the orientation (position angle, and for Blb-S, 
inclination), and assuming that the component of the magnetic 
field in the plane of the sk y lies between 19 fiG as derived by 
iMatthews & Wils^ d2002h and 100 ;uG, we constrained the an¬ 
gle p between each outflow and the magnetic fieldQto be/l > 36° 
for Bl b-N and/? > 29° for Blb- S. As shown in MHD simula¬ 
tions dCiardi & Hennebellell2010h . significant asymmetries are 


' The outflow orientation is defined by its projected axis posi¬ 
tion angle PA and inclination i. p can be obtained using cosifi) = 

B±cos(PA)cos{i)+B / / sin{i) 
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predicted between the red and blue lobes at such high inclina¬ 
tions, which may explain the shape of Blb-N and Blb-S out¬ 
flows. 


3.2. Physical properties of the outflows 


We extracted spectra of the four detected CH 3 OH lines at se¬ 
lected positions given in Table |2] and fitted these profiles with 
a two-component Gaussian, assuming the same velocity profile 
for all CH 3 OH lines. The narrow component centered near 6.5 
km s ' traces the quiescent component, while the broad compo¬ 
nent describes the molecular outflow. The line intensities were 
compared with predictions from the non-LTE radiative trans¬ 
fer c ode MADEX (ICernicharol20T^ using collision cross-sections 
from iRabli & Flowej ( 2010|) to derive n(H 2 ), Tk and CH 3 OH 
column densities N in the broad component. The uncertainties 
are 5 - 10 K for Tk, ~ 50% for ^(Hi), and ~ 20% for N, 
as shown by the scatter at nearby positions. Because the out¬ 
flows show structures down to the beam size in the integrated 
intensity maps (FiglTJ, we assumed that the length along the 
line of sight for the high-velocity gas is equal to the beam 
size, 2.3" or 530 AU, to derive a lower limit to the H 2 column 
density. This assumption led to CH 3 OH abundances relative to 
H 2 in the 10 ^ - 10 ^ range, which are intermediate between 
the high values, ~ 10 ~^, reached in chemically-rich outflows 
( iBachiller & Perez Gutierre z 1997h and the low values measured 
in quiescent gas. ~ 10^'° dGuzman et al.ll20l'^ . The H 2 densi¬ 
ties derived from MADEX range from 3 x 10^ to 9 x 10^ cm“^. 
They are comparable with the densities in the B ib envelope 
determ ined from the dust continuum emission by [Daniel et alJ 
(l2013h . Althoug h higher than the kinetic temperature of the qui¬ 
escent gas, 12 K (iFis et al.feOlOl) . the kinetic temperatures in the 
shocked gas remain moderate, between 20 and 70 K. 

Using the average properties adapted to each outflow, 
[CH 3 OH] = 7x10-* in Blb-S, and 10 * for Blb-N, we computed 
the mass in each lobe of the outflowing gas by integrating the 
emission in the velocity intervals listed in Tabled The resulting 
values are displayed in Table [3] together with other parameters. 
The outflow masses agree within a fac tor of two with the p revi- 
ous determinations based on CO(2-l) dHirano & Fiull20l'4i . ex¬ 
cept for the blue lobe of the Blb-S outflow, where the CO emis¬ 
sion is optically thick. The better sensitivity and angular resolu¬ 
tion allow us to obtain lower values for the dynamical time and 
outflow size, resulting in higher figures for the momentum and 
momentum flux by a factor of a; 5. The mechanical luminosities 
are at least a factor of ten lower than the bolometric luminosity 
of each source. 


3.3. Evolutionary stage of B1b-N and B1b 

Blb-N and Blb-S are located to the right of the class 0 YSOs 
in the diagram of normalized mom entum flux, Fc/L versu s nor¬ 
malized envelope mass MenvIL^'^ (iBontemps et al.lfl996h . with 
higher normalized envelope masses (> 1.1 and > 0.7) and a 
fairly similar normalized momentum flux (800 and 1600). Their 
luminosities are significantly lower than predicted by the scaling 
with either the envelope mass or the momentum flux, however; 
this is a typical feature of VeFFOs. B Ib-N and B I b-S lie in fairly 
massi ve cores with at least 0.36 Mq in 250 AU (iHirano & Fiul 
12014 . With such a mass reservoir, the final object is most likely 
a low-mass star. 

Previous theoretical studies have shown that the degree of 
turbulence, characterized by the ratio of kinetic to gravitational 


energy e, and the strength of the magnetic field, characterized by 
the ratio of the mass to flux ratio, to the critical mass to flux ratio, 
jU0, are two key parameters for the evolution of collapsing cores 
(e.g. Lloos et alluoT^ ICommercon et alJl2012^ . U sing the den¬ 
sity pr ofile and turbulent line width determined by [Daniel et al.l 
(I 2 OI 3 I) in Bib, we derive e ~ 0.3 and p - 2 -1 depending on 
the strength of the magnetic field. The physical parameters of 
the Blb-N outflow listed in Table [2 agree very well with those 
of the low outflow com ponent in numerical models of the earl y 
evolution of protostars (I.Toos et aT]l2012i 120131 : lMachidall20l4 . 
especially for highly magnetized cores. The significant degree 
of turbulence could e xplain the misalig nment with the magnetic 
field. As explained hv Lloos et al.l d2013 i. the localized rotational 
motions induced by turbulence lead to a preferred orientation 
for the collapsing core that is not related with the local direction 
of the magnetic field. It could also explain the different orien- 
tation of the Blb-N and B lb-S outflows. Models discussed by 
ICommercon et al.l (1201 2ah predicted the mass within 250 AU to 
be M 2502 it/ = 0.16 Mq for a total core mass of about 1 M© within 
4000 AU. These figures are compa tible with the total mass of 
Bib, ~ 2.5 Mq (Irianiel et al.ll2013l) and the mass of the com¬ 
pact mm sources, 0.36 M© each. The outflow masses, velocities, 
and mechanical luminosities predicted by the MHD models de¬ 
pend on p, with Mont - 3.3 x 10^* M©, Vma.x - 4.7 kms^', 
Lmech - 4.1 X 10^* F© for p - 2 . For p - 10, the out¬ 
flow mass increases to Mo„, = 2. x 10^^ M©, the luminosity to 
Lmech = 7 X 10^* F©, while the maximum velocity decreases to 
Vmax = 1.7 kms-'. 

These figures are computed about 1000 years after the for¬ 
mation of the FHSC. The agreement of the p = 2 case with 
the properties derived for Blb-N, including its short dynami¬ 
cal age, compact size, and cold SED, is excellent, supporting its 
early e volutionary stage. As already discussed by iHirano & Fiul 
(12014 . the more energetic outflow properties, higher luminosity, 
and warmer SED of Blb-S only poorly fit with the same mod¬ 
els, providing further support for the classification of Blb-S as a 
class 0 YSO, currently in a low accretion stage. 

All data so far are consistent with the classification of Blb- 
N as an FHSC, but do not provide definitive evidence. Higher 
angular resolution observations are necessary to probe the gas 
dynamics at the < lOOAU scale, and a more complete character¬ 
ization of the chemical and physical properties of B lb is required 
as well. 
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Appendix A: Data processing 

Appendix A. 1: Interferometric data 

We used the standard algorithms implemented in the software 
GILDAS/CLIC to calibrate the PdBI data. The radio-frequency 
bandpass was calibrated by observing the bright (12.8 Jy) quasar 
3C84. Phase and amplitude temporal variations where calibrated 
by fitting spline polynomials through regular measurements of 
two nearby (< 11°) quasars (3C84 and 03334-321). The PdBI 
secondary flux calibrator MWC 349 was observed once during 
every track, which allowed us to derive the flux scale of the in¬ 
terferometric data. The absolute flux accuracy is ~ 10%. 

To produce the continuum maps, we imaged and decon¬ 
volved the WIDEX data at 2 MHz resolution. This allowed us to 
identify all the detected lines and to remove them before synthe¬ 
sizing the continuum image. To subtract the continuum from the 
lines, we first synthesized continuum uv tables in a ~ 400 MHz 
frequency range devoid of lines close to the targeted line. This 
way, we did not need to take a potential variation of the contin¬ 
uum level with frequency into account. 

Appendix A.2: Singie-dish data 

Data reduction was carried out using the software 
GILDAS/CLAS^ . The data were fi rst calibrated to the T! 
scale using the chopper-wheel method dPenzias & BurrusI 19731) . 
The spectra were converted to main-beam temperatures (T^b) 
using the forward and main-beam efficiencies (Fetr and Befr) 
listed in Table 1. The resulting amplitude accuracy is 10%. A 
20 MHz-wide subset of the spectra was first extracted around 
each line rest frequency. We computed the experimental noise 

^ See http: //www. iram. f r/IRAMFR/ GILDAS for more information 
about the GILDAS softwares llPetvll200.'5l) . 


after subtracting a first-order baseline from every spectrum, 
excluding the velocity range from 4 to 9 km s ' LSR where the 
signal resides. A systematic comparison of this noise value with 
the theoretical noise computed from the system temperature, the 
integration time, and the channel width allowed us to Alter out 
outlier spectra (typically 3% of the data). The spectra were then 
gridded into a data cube through a convolution with a Gaussian 
kernel of FWHM~ 1 /3 of the IRAM-30m telescope beamwidth. 

Appendix A.3: Joint imaging and deconvoiution of the 

interferometric and singie-dish data 

Following iRodriguez-Fernandez et al.l (l2008h . the software 
GILDAS/MAPPING and the single-dish map from the IRAM-30m 
were used to create the short-spacing visibilities not sampled by 
the Plateau de Bure interferometer. In short, the maps were de¬ 
convolved from the IRAM-30m beam in the Fourier plane before 
multiplication by the PdBI primary beam in the image plane. Af¬ 
ter a last Fourier transform, pseudo-visibilities were sampled be¬ 
tween 0 and 15 m, which is the difference between the diameters 
of the IRAM-30m and the PdBI antennas. 

These visibilities were then merged with the interferometric 
observations. Each mosaic held was imaged and a dirty mosaic 
was built by combining these field s in the following optimal way 
in ter ms of signal-to-noise ratio jPetv & Rodriguez -Fernand^ 
I 2 OIOI) . The dirty intensity distribution was corrected for pri¬ 
mary beam attenuation, which induces a spatially inhomoge¬ 
neous noise level. In particular, noise strongly increases near 
the edges of the held of view. To limit this effect, both the pri¬ 
mary beams and the resulting dirty mosaics were truncated. The 
standard level of truncation was set at 20% of the maximum in 
GILDAS/MAPPING. The dirty image was deconvolved using the 
standard Hogbom CLEAN algorithm. The resulting data cube 
was then scaled from Jy/beam to the Tmb temperature scale us¬ 
ing the synthesized beam size (see Table 1). 

The H 2 CO emission covers most of the mosaic held of view. 
The emission structure thus seems to sit on a constant bright¬ 
ness that only depends on the frequency, not on the spatial po¬ 
sition. CLEAN deconvolution methods have many difficulties 
to properly deconvolve this “constant” emission. To avoid this, 
we subtracted the mean spectrum over this held of view of the 
single-dish data before processing them to produce the short- 
spacings. We then imaged and deconvolved the hybrid data set 
as explained above and added this mean spectrum back to the 
hybrid synthesis data cube (30m + PdBI) after deconvolution 
and conversion to the Tmb temperature scale. This treatment is 
correct because a constant emission is always resolved, that is, 
independent of the resolving power of the observatory. 

Appendix B: Spectral energy distribution 

We present in Pig. IB. II the spectral energy distribution of Blb- 
N and Blb-S and provide refined positions in Table IBT] The 
m easured continuum flu xes agree well with the v alues reported 
bv iPezzuto et al.l (1201 2l) and iHirano & Liul (l2014h . It is interest¬ 
ing to observe that Blb-S is more luminous than Blb-N at far- 
infrared and submillimeter wavelengths while the reverse is true 
long-ward of ~ 2.5 mm. The new PdBI data help to locate the 
crossing point of the spectral energy distributions. 
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Fig. B.l. Spectral energy dist ribution of Blb-S (blu e) and B Ib-N (red). 
The data points are taken from iHirano & Liul l l20l4) and include Spitzer, 
Herschel, and ground-based observations. The stars show the new PdBI 
measurements. 

Table B.l. Positions and fluxes of the protostellar sources 


Name 

RA 

Dec 

L 

(LG) 

Fi45 

(mJy) 

Blb-N 

03:33:21.21 

31:07:43.8 

0.15 

77 ±4 

Blb-S 

03:33:21.36 

31:07:26.4 

0.31 

96 ±5 

Blb-W 

03:33:20.30 

31:07:21.4 

0.17/0.25 

< 4 
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